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SUMMARY 


A method is given for determining the horizontal tail loads in 
maneuvering flight. The method is based upon the assignment of a load- 
factor variation with time and the determination of a minimum time to 
reach peak load factor. The tail load is separated into various com- 
ponents. Examination of these components indicated that one of the 
components was so small that it could be neglected for most conventional 
airplanes, thereby reducing to a minimum the number of aerodynamic 
parameters needed in this computation of tail loads. 

In order to illustrate the method, as well as to show the effect 
of the main variables, a number of examples are given. 

Some discussion is given regarding the determination of maximum 
tail loads, maximum pitching accelerations, and maximum pitching veloc- 
ities obtainable. 


INTRODUCTION 

* 

The subject of maneuvering tail loads has received considerable 
attention both experimentally and theoretically. Theoretically, methods 
and solutions have been derived for determining the horizontal tail load 
following either a prescribed elevator motion (references 1 to 3) or an 
assigned load— factor variation (reference V). 

The first approach has been adopted into some of the load require- 
ments where the type of elevator movement specified consists of linear 
segments whose magnitudes and rates of movement are governed by the 
assignment of a maximum initial elevator movement consistent with the 
pilot’s strength. The rates of movement and the time the elevator is 
held before reversing are so adjusted that the design load' factor will 
not be exceeded. 
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The results of reference 5 show, as is to he ejected, that only 
when the aerodynamic force coefficients are accurately known from wind- 
tunnel tests can good agreement he obtained between measured and 
calculated tail loads. At the design stage, however, only general 
aerodyna mi c and geometric quantities are available and some of the more 
important stability parameters are not known accurately. Thus, the work 
involved in the solution for the tail load following a given elevator 
motion is not considered to be in keeping with the accuracy of the results 
obtained. Consequently, there appears to be a need for an abbreviated 
design method of computing tail loads which, although incorporating 
approximations, will nevertheless be based on the theoretical considera- 
tions of the problem. 

If the load— factor variation with time is specified and the 
corresponding tail load, elevator angles, and load distributions are 
subsequently determined, a simpler and equally rational approach to the 
tail— load problem can be made. Although this approach has been used to 
a limited degree (reference 4), several shortcomings have limited its 
use. 

The purpose of this paper is to develop further the load— factor or 
inverse approach and to present a method of computing horizontal tail 
loads which is comprehensive and generally simple. To this end, (1) 
the shape of the load— factor curve and the minimrum time required to 
reach the peak load factor have been determined from an analysis of 
pull-up maneuvers that were available, (2) the minimum time required to 
reach the peak load factor has been determined from a theoretical analysis 
which is supported in some measure by statistical data obtained from a 
number of flight tests with airplanes of widely varying sizes, and (3) 
the equations relating the various quantities are presented. 


SYMBOLS 

b wing span, feet; also shape factor in equation (13) 

b-fc tail span, feet 

c chord, feet 

c mean aerodynamic wing chord, feet 

lift coefficient (L/qS) 

^m pitching-moment coefficient of airplane without 

horizontal tail (Mb/qS 2 ) 
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pitching-moment coefficient of isolated- horizontal- 
tail surface 

acceleration of gravity, feet per second per second 

2 

pitching moment of inertia, slug— feet 

radius of gyration about pitching axis, feet 

empirical constant denoting ratio of damping moment 
of complete airplane to damping moment of tail alone 

lift, pounds 

local lift at any spanwise station 

airplane mass, slugs (W/g) 

moment, foot-pounds 

airplane load factor at any instant 

maximum increment in load factor 

dynamic pressure, pounds per square foot 

wing area, square feet 
horizontal— tail area, square feet 
time, seconds 

time to reach peak of elevator deflection, seconds 
airplane true velocity, feet per second 
airplane weight, pounds' 

length from center of gravity of airplane to aerod ynam ic 
center of tail (positive for conventional airplanes), 
feet 

nondimens ional spanwise dimension 
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constants occurring in equations (13), ( 23 ) ^ (26), and 
(30) 

ponstants occurring in “basic differential equation 
(see equation ( 3 ) and table I) 

time to reach peak load factor, seconds 

mass density of air, slugs per cubic foot 

tail efficiency factor (q-fc/q) 

wing angle of attack, radians 

average angle of attack of horizontal stabilizer, 
radians 

tail angle of attack, radians 
angle of sideslip, degrees 
flight-path angle, radians 
attitude angle, radians (a + 7 ) 
elevator angle, radians 

downwash angle, radians 
tail setting, radians 



The notations d and 0, a and *0, and so forth, denote single 
and double differentiations with respect to t. 

Subscripts: 

o initial or selected value 


t tail 

max maximum value 

Z 0 zero lift 

geo geometric 


c 


camber 
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METHODS 

Method of Determining the Dynamic Tail Load 


Basic equations of motion .— The simple differential equations for 
the longitudinal motion of an airplane for any elevator deflection (see 
method given in reference 2) may he written as 


m.77 - 



qS 



Tj^qS^. AS = 0 



d ' C m a S 2 

da ^ h 


dC 


ALfeZt 


®fc 


dS 


s t 2 

V sr 


AS 


2 •* 

mky 9 = 0 


(2) 


Equations (l) and (2) represent summations of forces perpendicular 
to the relative wind and of moments about the center of gravity. (See 
fig. 1 for direction of positive quantities . ) Implicit in these equa- 
tions are the following assumptions: 

(1) In the interval between the start of the maneuver and the 
attainment of maximum, loads, the flight— path angle does not change 
materially; therefore, the change in load factor due to flight— path 
change is small. 

(2) At the Mach number for which computations are made, the aero- 
dynamic derivatives are linear with angle of attack and elevator angle. 

(3) The variation of speed during the maneuver may be neglected. 

(4) Unsteady lift effects may be neglected. 

• * • '* »« 

By use of the relations 0 = 7 + 0 ., 6 = 7+ a, and 6 = 7 + a, 
equations (1) and (2) are reducible to the equivalent second— order 
differential equation 


a + K^a + Kg Aa = AS ( 3 ) 

where K^, Kg, and are constants for a given set of conditions (see 

table I). 
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In equations (l) and (2), Aa, 7, 0, AS, and ALj. will, in a 

given maneuver, vary with time. Using the relations between 0, 7, a, 
and their derivatives permits equation (2 ) to he rewritten as follows to 
give the increment in tail load: 


^t = *T 


2 - 
inky cx. 


% 8 r,S. M 

x* as w-t 


In a still shorter form., equation (4) may he written as 


AT. . = AT . , + AL + + AL. + ALj. 

t '*a t c 


Equations (4) and (5) show that the tail— load increment (the increment 
above the steady— flight datum value) at any time is composed of four 

parts: AL, , associated with the angle— of— attack change; AL. , asso— 

% 

dated with angular acceleration about the flight path; AL^, associated 

with angular acceleration of the flight path; and AL^ , required to 

c 

compensate for the moment introduced by change in camber of the horizontal- 
tail surface. The load AL^ is generally small but in some extreme 

configurations may amount to 10 percent of the total increment and thus 
for the present it is retained in the development. 

If the load— factor— increment variation with time An is known, 
then by the usual definition 


so that 




^ w / s 
dCr 

a, = sJdM 

dCr 

— -** q 

da 

.. ii w/s 
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The following relation also exists between An and 7 : 

An g = jrY 


so that 



(7) 

(8) 


"When equations ( 6 ) to ( 8 ) are substituted into equations (^) and (5), 
the four tail— load components then become 


AL+ 

% 


dC. 


'm ¥S 


dC L**t 


An 


(9a) 


2 

= “ V 31 (9h) 

gSqx t ^r 


Wky 2 . 

n (9°) 

dC_ g 2 

^0 * vr ’*«■ ^ (9 4) 


.Thus, if the variation of the load factor with time An and the 
geometric and aerodynamic characteristics of the airplane were known, 
the first three components of the tail load could be found immediately. 

The magnitude of the fourth component, that due to horizontal— tail camber, 
would follow from equation (3) in which the elevator angle is seen to be 


A 6 = g- + 

E 3 K 3 


. *2 
a + =- Ax 

K 3 


( 10 ) 


Substitution into equation (10) of the values of Ax, d, and a, 
from equation ( 6 ) yields the value of the elevator angle at any instant 
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A6 = 


W/S 


dO- 

£ 3 5T’ 


^n + K^h + K 2 An) 


(ID 


so that, finally, the fourth component is given as 


AL , = T). -- ^ S — (n + E^h + K 2 An) 

t c <3S 't b^ d^ V 1 ^ f 


( 12 ) 


T - K, 

da 3 


The procedure outlined shows that the tail— load magnitude can he 
determined if the load— factor variation is known. 

Types of load— factor variation .— The relation between the tail load, 
the geometric and aerodynamic characteristics, and the load factor having 
been established, it is desirable to establish a load— factor variation 
which is reasonable as well as critical insofar as loads are concerned. 

The maximum value of load factor is usually specified; however, there 
are many possible variations for the shape. Begardless of the details 
of shape, the load factor may be considered to rise smoothly and con- 
tinuously to a maximum, the rate of rise depending upon several variables. 
Beyond the maximum value of the load factor the return to initial condi- 
tions can, at the will of the pilot, be either gradual or rapid. 

Experiments as well as theoretical studies have already indicated 
that the maneuver that combines maximum angular and linear accelerations 
causes critical loads in both the wing and tail. One such maneuver occurs 
when the maximum load factor is reached as rapidly as possible by using 
an initial elevator movement which is greater than that required to reach 
a given steady— trim value of the load factor. Tills initial elevator move- 
ment is followed by a rapid checking of the maneuver either by returning 
the elevator quickly to neutral or by reversing the controls. 

The shape of the load— factor curve for such a maneuver may be 
expressed approximately by several analytic functions, one of which is 

An = at^e -0 ^ (13) 


By way of illustration, figure 2 shows details of the shape of the load- 
factor curve obtained with the use of equation ( 13 ) for which the constants 
have been adjusted so that an 8g peak is reached in 1 second. By further 
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adjustment of the constants the load factor can, within certain limits, 
be made to rise to any specified peak and to diminish in any prescribed 
manner. 

Because the posit ire slopes obtained from equation (13) are always 
greater than the negative slopes, the positive angular accelerations are 
greater than the negative ones. In general, this condition is true for 
most high g critical maneuvers performed by most classes of airplanes, 
but maneuvers may occasionally be performed for which the reverse may be 
true, particularly for small airplanes. 

Determi nation of constants. — i'rom equations (9), (11), and (12) 
the required quantities relating to load factor are seen to be An, n, 
and ii. Since the increment An is to be given by 


An = at b e“ ct 


then at maximum load factor 



Thus t = £ 
c 


so that 



at maximum load factor. Let N = A n^y . Then 




(13) 


(14) 


(15) 


(16) 


(17) 


Equation ( X7) is in nondimens ional form where X is the time to reach 
the peak load factor and b is a constant. 
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When equation (17) is differ satiated., the first and second deriva- 
tives become 

W> 

and 

hA. £ _ £n ^2 
H N 

In equations (17) to (19) the quantities IT, A, and b are now 
required in order to determine the variation of An, h, and n. The 
value of IT is immediately available from the required maneuver load 
factor, whereas the time to reach the peak load factor X can be ob- 
tained from examination of available records or by specification. The 
constant b, as may be seen from equation (17), can best be described 
as a "shape" factor and has no particular physical significance. 

The values of X and b should be associated with a maneuver which 
produces maximum tail loads. Therefore the time X to reach peak load 
factor should be the minimum possible consistent with possible pilot 
action and airplane response. The shape factor b should also be con- 
sistent with both of these. 

In connection with the determination of the minimum time to reach 
peak load factor, the results shown in figure 3 for a typical airplane 
are informative. Figure 3(a) shows the load— factor variation following 
several abrupt jump elevator movements. The load factor varies with the 
elevator position, but the time to reach peak load factor does not. 

Figure 3(t) shows the load— factor variation for several abrupt hat- 
shape elevator impulses. Again the load factor is seen to vary with 
the amount of elevator deflection but the time to reach the peak value 
remains constant. Although the time to reach the peak load factor 
shown in figure 3(t) remains constant, it is seen to be less than that 
shown in the previous .case; therefore, an impulse elevator motion pro^ 
duces a smaller value of A, than the jump type. 

Because of inertia and elasticity in the control system, the pilot 
cannot move the elevator instantaneously but requires some finite time 
tj_ to do so. A possible critical type of elevator impulse thus appears 

to be one which increases linearly to maximum and decreases at the same 
rate to zero. In order to determine the minimum time, to reach peak load 
factor associated with such a variation, the equation of motion (equa- 
tion (3)) has been solved for the triangular elevator impulse for 
airplanes of various static stabilities and damping. 
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The results of the computations are given in figure A in which the 
minimum time X to reach peak load factor is plotted against the time t 
required to deflect the elevator. 

For completeness the curves of figure it- are labeled for the actual 
values of Kg employed in the confutation as well as for relative values 
of stability. By a series of computations the damping term Kp was 
found, as was to be expected, to have only a secondary effect on X- 
The curves apply to an average value of the damping constant. The upper 
curve, labeled "low stability,” should be associated with rearward 
center— of-gravity positions (that is, low static margin) in combination 
with one or both of the following: low dynamic pressure or heavy air- 

planes. The lower curve, labeled "high stability," would be associated 
with forward center— of-gravity positions in combination with one or both 
of the following: high dynamic pressure or light airplanes. It is seen 

that X increases almost linearly with tp and also increases when the 
restoring forces are reduced, that is, when the stability is reduced. 

A preliminary value of the shape factor b (required in equa- 
tions (17) to (19)) was initially determined from flight records of 
typical impulse maneuvers by fitting curves of the type given by equa- 
tion (13) through several points of the actual time histories and 
determining the constants. The results of this first step were then 
modified by the results of the same computations which had been made to 
determine X, and the variation of b with tp given in figure 5 was 
obtained. Since the b factor is not found to be critical, an average 
value of 5.0 is suggested, although as a refinement the values from 
figure 5 may be used. 

The question of the value of tp to use is one which must be 
solved either from experience or from a knowledge of the characteristics 
of the controls and the control system. For conventional airplanes 
having the usual amounts of boost and no rate restrictors, the following 
values of tp are suggested as representative: 


Fighters or small civil airplanes with weight limit from 

about 500 to 12,000 pounds , seconds 0.20 

Two— engine airplanes with weight limit from 25,000 

to 45,000 pounds, seconds 0.25 

Four-engine airplanes with weight limit from 50,000 

to 80,000 pounds, seconds . . ........... 0.30 

Airplanes with weight limit above 100,000 pounds, seconds . . ; 0.40 
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Tlie Minimum time X given in figure 4 was actually established 
.separately from the adopted load— factor variation; therefore, in applying 
the inverse method, the derived elevator impulse would not he expected 
to agree in detail with the "tent" type impulse used in the derivation. 

The first three tail— load components can now he computed hy insert- 
ing the values of An, h, and n from equations (17) to ( 19 ) into 
equation ( 9 ) and using appropriate values of X from figure 4 . In 
order to facilitate this computation, curves of An/u, hx/H, and nX 2 /lf 
plotted against t/x are given in figure 6 for the suggested value of 
h = 5. Actually to apply the results of figure 6 it is convenient to 
find first the components AL, , and so forth, in terms of the nondi— 

a 

mensional time t/x and then to convert to time t in seconds. In 
order either to compute the fourth component or to obtain the elevator 
angles for use in chord loading, the constants Kj_, 'Eg, and E3 of 
equation (3) must also he known. 


Thus, in terms of t/x and the ordinates of figure 6, the various 
tail— load components are 


AL+ = —S’ ^^(ordimate of fig. 6(a)) 
dC^ hx t \ > 


(20a) 


-W 2 ^ 2 H / 

dC, *2 V 
eSq*t“a£ 


Ordinate of 


fig. 6(c)) 


- Wkv 2 
AT. 1- = - 

7 


■2— ^Ordinate of fig. 6(h)) 


(20b) 


(20c) 


AL+ = 


dCm fe s t 2 W 
dS t ^t^t S 


H 

Ordinate of fig. 6(c) 

dC T 

Kr-sS 

[ x 2 


K]_ ^Ordinate of fig. 6(h)) 
" ' ~X 


+ Kg (ordinate of fig. 6(a)) 


( 20d) 


The constants K^, Kg, and defined in table I are the same as those 

given in reference 2, except for changed signs caused hy specifying x-fc 
as positive. 
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/ 

The conversion to time t is made by multiplying values of the 
base scale t/l by X. 

Sample calculations for incremental tail loads .— The results of 
several examples are given to illustrate not only the method but also 
the effect of each of a number of variables on the incremental tail load 
of a typical fighter airplane, the geometric and aerodynamic character- 
istics of which are given in the following tables. In order to illustrate 
the effect of static stability, results have been computed for three 
center-of— gravity positions with the assumption that an 8g recovery is 
made at 19,100 feet from a vertical dive at an equivalent airspeed of 
400 miles per hour. In order to illustrate the effect of the time of 
the elevator impulse on the tail load, computations were carried out at 
one of the center— of -gravity positions for several values of t^. The 

cases considered and the airplane characteristics follow: 


GEOMETRIC CHARACTERISTICS 


Gross wing area, S, square feet . . . . . . . . . 300 

Gross horizontal— tail area, S^, square feet .......... 60 

Airplane weight, W, pounds ................... 12,000 

Wing span, b, feet ............ « . 4l 

Tail span, b^, feet 1 6 

Radius of gyration, ky, feet 6.4 

Distance from aerodynamic center of airplane less tail to 
aerodynamic center of tail, x^., feet: 

Center of gravity, 29 percent M.A.C. ............ 20.0 

Center of gravity, 24 percent M.A.C. ............ 20.3 

Center of gravity at aerodynamic center . . ........ 21.0 


AERODYNAMIC CHARACTERISTICS 


Slope of airplane lift curve, dC^/dcc, radians 
Slope of tail lift curve, dCj^/da^, radians 

Downwash factor, de/da . . . 

Tail efficiency factor ^qyq 

Empirical airplane damping factor, K 

Elevator effectiveness factor, radians 

Rate of change of tail moment with camber due to 
elevator angle, dC^/dS, radian ....... 



4.87 

3.15 

0.54 

1.00 

1.1 

1.89 


-0.57 
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Rate of change of moment coefficient with angle of attack 


for airplane less tail, dC m /da, radians: 

Center of gravity, 29 percent M.A.C. ........... 0.625 

Center of gravity, 2k percent M.A.C. ....... . . . . 0.1*03 

Center of gravity at aerodynamic center .......... 0.000 


The specified conditions for the sample computations are given in 
table II. The computed results for tail components are given in fig- 
ures 7 and 8. Figure 7 gives results for varying the center of gravity 
and figure 8 gives similar results for varying tp. The tail— load 
components are computed from equations (20) and the derived elevator 
angles from equation (ll). If the increment in tail load due to camber 
and the incremental elevator angle are not required, the K values need 
not be computed and the computations are considerably shortened. Fig- 
ures 7 and 8 show that a maximum error of only about 4 percent is 
introduced by this omission.- 


Method of Determining the Total Tail Load 


The initial or steady— flight tail load and elevator angles to which 
the computed incremental values are to be added must also be determined. 
In steady flight, the horizontal tail furnishes the moment required to 
balance the moments from all other parts of the airplane so that the 
initial load may be written as 


CVqSc 

- - V 


dC m WS 
3% bx t 


cos 


7o 


( 21 ) 


Thus the total tail load at any time in a maneuver is composed of the 
four previously mentioned parts plus the components given in equa- 
tion (21). Only the first term of equation (21) represents a new type 
of load because the second term is a load of the type given by equa- 
tion (9a) or equation (20a) and its effect may be immediately included 
in the computations by multiplying the ordinates of figure 6 by 
R + cos 7 q instead of by R. 

The initial elevator angle required to balance the airplane in 
steady flight varies with airplane Cp, and center-of— gravity position 
so that, in general, 5 0 must be obtained from wind-tunnel data. With- 
out results of wind-tunnel tests, a rough rule which can be used as a 
guide at the design stage in determining the elevator position is that 
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the final elevator setting will he so adjusted by repositioning of the 
stabilizer setting during acceptance tests that it will be near a zero 
position at the cruising speed and at the most prevalent center— of— 
gravity position. 


Method of Determining Maximum Yalues 

Maximum tail loads and angular accelerations .— The method outlined, 
enables a point— by-point evaluation to be made of the quantities that 
determine the tail load. Such detail may often be unnecessary and the 
procedure may be shortened by evaluating only those points near the load 
peaks or, alternatively, by accepting an approximation to the results. 

One such approximation which may be made is to balance the airplane at 
the combinations of load factor and angular acceleration which would 
result in maximum up and down tail loads. 

Figure 7 shows that the maximum down tail load in a pull— up occurs 
near the start of the maneuver and before appreciable load factor is 
reached. This maximum load is practically coincident with the negative 
maximum in the Lx. tail— load component. Since, for a given configu— 

a 

ration, this component increases as the center of gravity is moved forward 
and since the steady— flight down load increases with speed, the maximum 
down tail load in a pull-up occurs at the highest design speed in combin- 
ation with the most forward center— of-gravity position. 

Figures 7 and 8 show that at the time of the maximum down-tail— load 
increment the elevator is near but has not quite reached its peak position. 
Also at the time of maximum up-tail— load increment the elevator is near 
its zero position, although it may be on either side of this position 
depending upon the stability and the time t-j_. These results suggest 
that the maximum down load for the elevator and the hinge brackets would 
occur with the airplane center of gravity well forward and at the start 
of the maneuver . The maximum load for the stabilizer is likely to occur 
at the peak load factor. 

Figure 7 also shows that the up tail load occurs near the peak of 
the L-fcj component as well as near the positive maximum peak in the 

L-fc„ component. Since the component increases as the center of 

gravity is moved rearward and since a decrease in speed generally reduces 
the initial down load, the maximum up tail load occurs at the upper left- 
hand corner of the Y— n diagram for the most rearward center— of-gravity 
position. 
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The maximum tail load in a pull— up maneuver may he written as 


“max 


VSC + Si (n + cos 7o ) 


x t dC L hx t 


10 

*t 


( 22 ) 


where the sum of the second and third terms is to he a maximum in the 
maneuver. From the previous discussion the load— factor increment at 
maximum down load is nearly zero and at maximum up load it is^nearly 
equal to N so that if the positive and negative values of 0 ^ can 

he determined, a relatively simple method for determining maximum, loads 
is available. 

Since hy definition 0 = a + 7 , an expression for angular acceler- 
ation can he derived from equations ( 6 ) and ( 7 ) and written in the form 


-- W/S •• 


a 



The maximum angular acceleration can he approximated hy 


v = B .2 IS- a 

°max n ^2 

da ^ 


+ 


° YX 


(23) 


For the maximum positive pitching acceleration, B is the maximum 
positive ordinate in figure 6 (c) and C is the ordinate of figure 6 (h) 
at a value of t/x for which B was determined. Thus, B is 6.5 and 

C is 0-95 for this example. 

For the maximum negative pitching acceleration, B is the maximum 
negative ordinate in figure 6 (c) and C is the ordinate of figure 6 (h) 
at a value of t/x for which B was determined. Thus, B is —5-8 and 

C is 0.80. For use in equation ( 23 ) the values of X for the maneuver 

are available from figure 4 and the other quantities are available from 
the conditions of the problem. The maximum loads can he given hy the 
following equations: 
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For maximum up tail load in the pull-up: 


L+ ■'* — + ^ (K + 1.0) + — f fe, 8 

t max+ *t dCt bx t v x t X dC^ v 


l da 


q.X 


(24a) 


For maximum down tail load in the pull-up: 


= ¥ + + < £to > 

\3T / 


For push-downs to limit load factor, equations (2l}.a) and (21b) still 
apply with changed signs for N and changed directions- for and 

max+ 

L+ . A question arises as to whether the ma/H mum down tail load at 
max— 

the start of a pull-up with forward center— of-gravity position is greater 
than that which would occur when pulling up from a negative load— factor 
condition, with the center of gravity in the most rearward position. 

This can be determined only by computing both cases and seeing which is 
the larger. 


Maximum value of angular velocity .— The maximum value of the pitching 
angular velocity in the pull-up may also be found in a manner similar to 
that used to obtain the maximum angular acceleration. Since 6 = a + y 
and the relations involving these quantities in terms of load factor are 
given by equations (6) and (7), the following equation may be written: 


_ n W/S An g 
7 

da y - 


( 25 ) 


The maximum angular velocity may be approximated by 
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where D is the maximum positive ordinate in. figure 6(h) and. E is the 
ordinate of figure 6(a) at a value of t/X for which D was determined.. 
Thus D, for this example, is 1.95 and E is 0.48. 

In the steady turn or pull-up at constant g, the angular velocity 
is usually given by the expression 0 = 1.0 The difference between 

the factor 1.0 of this expression and the factor 0.48 of equation (26) 
is more than made up by the angle-of— attack component of the angular 
velocity. 


Approximate Method of Determining Load Distribution 


Symmetrical loading .— The spanwise distribution of the total load 
can be formulated with various degrees of exactness. If information 
regarding details of the angle— of-attack distribution across the span 
were available, then an exact solution could be obtained for the loading 
with the use of existing lifting— surface methods. The following method 
may be used as a first approximation to the solution. 

Erom the total tail load, the total tail lift coefficient Gt can 

readily be found. The average effective angle of attack a of the 
stabilizer portion is given in the definition 

\ i iy * + Jo + AS ) § (27) 


where only a is assumed as unknown and c-t and c 7 may be taken 

as the rates of change of section lift coefficient with a and 5, 
respectively. 

Thus, for constant elevator angle across the span. 



(28) 
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In a practical case 150111 integrals in equation (28) need be evaluated 
only once for a given configuration and Mach number. A plot of d 

against with 6 as a parameter would he useful in further compu- 

t 

tations. With a known as a function of Cr and S, the local lift 

L t 

at any spanwise station is then obtained from the expression 


l = c 


t^c = [c 2a a + c 2& (S 0 + AS)] qc 


(29) 


Uhs ymmetr i cal loading .— % to this point the total loads have been 
assumed to be symmetrical about the airplane center line, whereas, in 
reality, the load may have an unsymmetrical part. The sources of this 
dissymmetry may be due to uneven rigging, differences in elasticity 
between the two sides, or to effects of slipstream, rolling, and sideslip. 
The first two sources are usually inadvertent ones while the last two are 
difficult to determine without either wind-tunnel tests or a knowledge of 
how the airplane will be operated. Present design rules regard in g dis- 
symmetry of tail load are concerned more with providing adequate design 
conditions for the reqr of the fuselage than with recognizing that at 
the maximum critical tail load some dissymmetry may exist. 

Tests in the Langley full-scale tunnel (reference 6) and flight 
tests (reference 7) of a fighter— type airplane, as well as unpublished 
flight tests of another fighter-type airplane, indicate that the tail- 
load dissymmetry varies linearly with angle of sideslip so that the 
difference in lift coefficient between the two sides of the tail can be 
given as 


C L — Ct = A3 ' (30) 

jRight Uieft 

The average values of A per degree found for the two fighter— type 
airplanes are approximately 0.01. Wo similar values are available for 
larger airplanes nor for tail surfaces having appreciable dihedral. 

In maneuvers of the type considered herein it is doubtful that 
angles of sideslip larger than 3° would be developed at the time the 
maximum tail load is reached. If the value of the sideslip angle at 
the time of maximum tail load can be established, equations (27) to (29) 
are easily modified to include this effect, provided the approximate 
value of A is known. 
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Chordwise loading . - The chordwise distribution can he determined 
for any one spanwise station in either of two ways. One way for design 
work is outlined in reference 8. A knowledge of the airfoil section 
and the quantities contained in equation (29) suffices for this deter- 
mination. 

If pressure-distribution data are available for a similar section 
with flaps, an alternate way would, he to distribute the load chordwise 
according to the two-dimensional pressure diagrams with the use of the 
computed values of section lift coefficient and elevator angle. 


DISCUSSION 


' The method presented is another approach to the determination of 
tail loads, Ih-om the results given in figures 7 and- 8, it can he seen 
that the camber component is so small that for all practical cases 

it may he omitted with considerable simplification in the computation of 
tail loads. This omission reduces to a minimum the number of aerodynamic 
parameters needed to compute the tail loads. 

It is possible, in the application of the present method with the 
use of the suggested values of tj_, that the derived elevator angles 

may not be within the pilot's capabilities. Since it must be assumed 
that all airplanes, to be satisfactory,, should have sufficient control 
to reach their design load boundaries, such an occurrence requires only 
that the time to reach elevator peak deflection tq be increased so as 

to reduce the elevator angle. The results of figure 8, in which the 
time tq is varied, furnish a useful guide for determining the in- 
crease tq that might be required. 

If sufficient information is available, it is recommended that 
existing lifting— surf ace methods be used in determining the spanwise 
distribution of the total load; however, if Information of the angle- 
of-attack distribution across the span is not known, the method presented 
may be used as a first approximation. 

Along some of the boundaries of the V-n diagram, tail buffeting may 
occur. Measurements show that buffeting usually occurs along the line 
of maximum lift coefficient and again along a high-speed buffet line 
which is associated with a compressibility or force break on some major 
part of the aiiplane, All aiiplanes are subject to buffeting at the 
design conditions associated with the left-hand comer of the V— n 
diagram. 
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Only high-speed and/or high-altitude airplanes are capable of reaching 
the other boundary. Measurements show that the oscillatory buffeting 
loads may be so high that the designer should at least be cognizant of 
them at the. design stage. 

The maximum angular acceleration varies inversely with airspeed 
and directly with the load factor, with the contribution due to acceler- 
ation in angle of attack likely to be more important than the angular 
acceleration of the flight path. A somewhat similar variation is indi- 
cated for the maximum angular velocity (equation (26)) where it is seen 
by direct substitution that the part due to angle of attack is likely 
to be larger than the part due to the angular velocity of the flight 
path. 


CONCLUDING KMARE3 


A simple method has been presented for determining the horizontal 
tail loads in maneuvering flight with the use of a prescribed incremental 
load— factor variation. 

The incremental tail load was separated into four components repre- 
senting a, a, y , and c. The camber component L^. is so small that 

c 

for most conventional airplanes it may be neglected, thereby reducing to 
a minimum the number of aerodynamic parameters needed in this computation 
of tail loads. 

An approximate method is presented for predicting maximum angular 
accelerations and maximum angular velocities . 

The method indicates that maximum tail loads in a pull-up occur at 
forward center— of-gravity positions and early in the maneuver. The 
maximum down tail loads in a pull-up occur at the highest design speed 
in combination with the most forward centeixif— gravity position. The 
maximum up tail load occurs at the upper left-hand corner of the V— n 
diagram for the most rearward center-of— gravity positions. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va. , February 9, 1950 
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TABLE I 
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TABLE U 


SPECIFIED COHDITIONS OE SAMPLE PROBLEM 


Increment in load factor . . . . . . . . • • . • . . . « . 8.0 

Altitude, feet . . . . . * • • • . 19,100 

Air density, sing per cubic foot . , .. . ....... . . . 0.001306 


Case 

Gg 

(percent M.A. C. ) 

t l 

% 

Kg 

e 3 

X 

(fig. !<-) 

1 

£L • C * 

0.2 

4.93 

30.4 

-33.^ 

0.45 

2 

2k 

.2 

k.72 

16.2 

-32.2 

.50 

3 

29 

.2 

i(-. 6 l 

8.45 

-31. 7 


k 

2k 

,k 

k.12 

16.2 

-32.2 

■ 77 

5 

2k 

.6 

k-72 

16.2 

-32.2 

1.02 
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Figure 1.- Sign conventions employed. Positive directions shown. 
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0 Z 4 .6 .8 1.0 IZ 1.4 1.6 /.8 

Time , vsec 

Figure 2.- Variation of load-factor increment. An = N250t5*53 e -5>53. 
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Figure 5*- Variation of shape factor b with t 
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Figure 6.- Variation of incremental-load-factor curves with time ratio. 
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